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RESEARCH  OBJECTIVES 


The  objective  of  the  research  is  to  attempt  to  develop  systematic 
techniques  for  determining  the  response  to  dynamic  impact,  when  the  shape  of 
the  impacting  object  is  known.  It  is  hoped  that  the  dynamic  response  can  be 
obtained  by  the  use  of  techniques  that  accurately  model  the  stresses  near  the 
point  of  impact  as  well  as  the  stresses  away  from  the  point  of  impact.  This 
technique  is  one  that  involves  a  global-local  technique,  where  a  more  exact 
description  of  the  stress  distribution  is  made  near  the  contact  (for  example, 
elasticity  theory)  and  a  more  approximate  theory  is  used  away  from  the  point 
of  load  application  (beam,  plate  or  shell  theory).  An  analysis  should 
incorporate  tne  following  features: 

a.  The  elastic  and  or  plastic  stresses  in  the  vicinity  of  the  impacting 
object  must  be  »,;c  jrately  cnaracteri zed . 

o.  An  accjrjt^  r  c '  m  or  tne  material  must  be  used,  especially 

wren  pons  ’  Ce~ :  n  je  near  the  impact. 

c.  At  points  distant  ‘>on  tne  load  the  analysis  must  accurately  model 
tne  dynamic  response.  Effects  of  support  members,  joints  or  other 
restraints  on  tne  response  must  be  carefully  modeled. 

d.  An  evaluation  of  the  potential  towards  damage,  such  as  cracking, 
plastic  deformation,  delamination,  etc.  must  be  made. 

The  research  objectives  include  the  study  of  both  static  and  dynamic 
solutions  to  indentation  problems.  The  static  ones  are  solved  in  order  to 
understand  the  elasticity  result  that  will  tend  to  form  the  basis  for  the 
local  solution.  This  solution  will  be  incorporated  into  a  dynamic  structural 
mechanics  solution  for  a  beam,  plate  or  shell  to  obtain  results  of  both  a 
global  and  local  character.  The  first  year  of  research  concentrated  on  the 
solutions  appropriate  to  beam  theory  and  to  axially  symmetric  plate  theory. 


STATUS  OF  THE  RESEARCH 


The  research  Is  progressing  along  the  lines  of  the  original  proposal. 

The  following  are  problems  that  involve  the  global-local  approach  that  have 
been  solved  for  the  static  case: 

1.  SMOOTH  INDENTATION  OF  AN  ISOTROPIC  CANTILEVER  BEAM  by  L.M.  Keer  and 

W.  Schonberg.  The  elastic  response  of  a  cantilever  beam  of  finite  length 
to  a  frictionless  cylindrical  and  flat  indenter  is  studied.  Solutions 
are  obtained  through  a  global-local  technique,  which  accounts  for  both 
the  local  behavior  near  the  indenter,  as  well  as  the  global  beam 
behavior.  Results  are  compared  to  Hertz  theory  and  to  beam  theory 
solutions  and  comments  on  the  validity  of  the  linear  elastic  analysis  are 
made . 


2.  SMOOTH  INDENTATION  OF  A  TRANSVERSELY  ISOTROPIC  CANTILEVER  BEAM  by  L.M. 
Keer  and  W.  Schonberg.  The  plane  elastic  response  of  a  transversely 
isotropic  beam  of  finite  length  to  a  frictionless  cylindrical  and  flat 
indenter  is  studied.  Local  indenter  stresses,  as  well  as  displacements 
and  rotations  are  computed  for  each  case  and  plotted  for  various  ratios 
of  contact  width  to  beam  length,  and  for  various  positions  of  the 
indenter.  Each  problem  is  solved  using  a  nearly  isotropic  and  a  highly 
anisotropic  material.  Results  are  compared  to  those  of  the  isotropic 
study,  and  to  Hertz  theory  and  beam  theory  solutions. 

An  extension  of  these  problems  to  the  dynamic  case  is  presently  unden 
way.  Possible  titles  are  given  in  "PAPERS  IN  PREDARATION" 

The  following  are  toe  dynamic  problems  that  have  been  solved  and  are  presently 
in  press  or  being  considerec  for  publication. 
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elastic  plate  was  solved.  The  indenter  was  assumed  to  have  a  radius  of 
curvature  sufficiently  large  that  areas  of  contact  of  the  order  of  the 
plate  thickness  could  be  considered.  Clamped  and  simply  supported  cases 
of  plate  impact  were  considered. 

2.  LOW  VELOCITY  IMPACT  ON  A  TRANSVERSELY  ISOTROPIC  PLATE  by  L.M.  Keer  and 

K.T.  Woo.  The  problem  of  the  dynamic  impact  of  a  rigid  indenter  striking 
a  transversely  isotropic  circular  elastic  plate  was  solved.  The  problem 
geometry  and  formulation  was  similar  to  (1)  except  that  the  material 
properties  were  for  transverse  isotropy.  Three  materials  were  chosen  for 
this  case:  cadmium,  magnesium,  and  a  laminated  composite. 

In  addition  to  the  above  papers  the  dynamic  impact  of  initially  stressed 
beams  and  cantilever  beams  are  papers  in  preparation.  They  should  be 
submitted  for  publication  relatively  soon. 
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SMOOTH  INDENTATION  OF  AN  ISOTROPIC  CANTILEVER  REAM 


by 

Leon  M.  Keer  and  William  P.  Schonberg 
Department  of  Civil  Engineering 
Northwestern  University 
Evanston,  IL  60201 


Abstract 


The  response  of  an  isotropic  cantilever  beam  of  finite  length  under  the 
action  of  frictionless  cylindrical  and  flat  indenters  is  studied.  Solutions 
are  obtained  through  a  local-global  technique,  which  accounts  for  both  the 
local  behavior  near  the  indenter,  as  well  as  the  global  beam  behavior.  The 
method  of  analysis  superposes  an  infinite  layer  solution  derived  through  the 
use  of  integral  transforms  with  a  pure  bending  beam  theory  solution.  Local 
indenter  stresses,  as  well  as  displacements  and  rotations  are  computed  for 
each  case  and  plotted  for  various  ratios  of  contact  width  to  beam  length,  and 
for  various  positions  of  the  indenter.  Where  possible,  the  results  are 
compared  to  Hertz  theory  of  contact  stresses  and  to  beam  displacement  and 
rotation  solutions. 


1.  INTRODUCTION 


Consider  the  finite  layer  of  length  l  and  thickness  k  shown  in  Fig.  1. 

The  layer  is  fixed  at  one  end,  free  on  the  other,  and  loaded  by  an  indenter 
centered  at  a  distance  £q  from  the  fixed  end.  Although  the  theory  is 
applicable  for  an  arbitrary  geometry,  the  types  of  indenters  studied  are 
cylindrical  (Fig.  la)  and  flat  (Fig.  lb).  Such  problems  are  seen  quite  often 
in  mechanical  applications  and  may  also  serve  as  models  for  impact  phenomena 
in  gears  and  turbine  blades.  Usual  methods  which  use  a  beam  theory  solution 
to  obtain  an  overall  load-displacement  relationship  and  then  a  Hertzian 
contact  solution  to  calculate  local  stresses  under  the  indenter  are  fairly 
limited,  as  'will  be  shown  in  this  paper.  The  method  of  solution  used  in  this 
paper  is  the  superposition  of  an  elasticity  solution  with  a  beam  theory 
solution.  In  this  way,  local  contact  stresses  are  represented  by  the 
elasticity  solution,  while  the  global  stresses  are  represented  by  the  beam 
theory  solution. 

A  problem  of  the  type  considered  here  has  already  been  solved  by  Keer  and 
Miller  flj,  and  by  Keer  and  Ballarini  [2]  for  a  beam  that  is  clamped  or  simply 
supported  at  both  ends  and  loaded  by  a  cylindrical  indenter.  Although  not 
treated  numerically  here,  the  intermediate  case  of  a  partially  fixed  end  can 
be  solved  for  each  type  of  indenter  in  a  similar  manner.  Furthermore,  the 
methods  used  in  this  paper  can  be  modified  and  expanded  to  study  plate  or 
shell  oroblems  involving  similar  loadings  [31. 

The  physical  quantities  of  interest  are  the  stress  distribution  under 
°ach  indenter,  and  the  displacement  and  rotation  of  the  beam  under  each  of  the 
i  ^or  ‘e"S .  ~he  ratio  of  contact  wi  nth  to  beam  thickness,  c'h,  is  t  reaped  an  a 
krq-wn  pa-ar-eter  •, n  hot n  probl ems .  I-  the  cylindrical  i ".’enter  ere  hi o-n, 
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ratio  of  indenter  arm  length  to  beam  thickness,  ?.g/h,  i s  treated  as  an 
additional  known  parameter.  Mowever,  in  the  flat  indenter  problem,  the 
parameter  Ai/h  =  £0/h  -  c/h  (i.e.  the  distance  from  the  fixed  end  of  the  beam 
to  the  left-hand  end  of  the  indenter)  is  introduced  as  an  additional  known 
parameter.  In  this  manner,  for  a  given  £i/h,  when  the  quantity  c/h  is  varied 
over  a  range  of  values  the  phenomenon  known  as  "receding  contact"  can  be 
observed  in  the  flat  indenter  problem.  Ry  varying  the  proper  parameters  as 
required  by  the  type  of  problem,  then  calculating  deflections,  rotations, 
loads  and  stresses,  an  extensive  study  of  the  response  of  an  isotropic 
cantilever  beam  subjected  to  various  loading  conditions  can  be  performed.  The 
results  obtained  are  compared  to  beam  theory  solutions  and  Hertz  contact 
solutions  for  accuracy  assessments. 


2C  CASE  I:  CYLINDRICAL  INDENTER  -  PROBLEM  FORMULATION 


The  problem  to  be  solved  is  that  of  an  elastic  layer  of  thickness  h  and 
length  £  indented  by  a  cylindrical  punch  on  its  upper  surface  (Fig.  la).  The 
conditions  at  the  ends  of  the  layer  are  those  of  a  cantilever  beam,  clamped  on 
me  left  and  free  on  the  right.  The  solution  of  the  problem  will  be  achieved 
by  a  suitable  superposition  and  matching  of  an  elasticity  solution  and  a  beam 
theory  solution. 

The  boundary  conditions  for  the  elasticity  problem  can  be  written  as 
fol 1 ows : 
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The  beam  theory  boundary  conditions  for  the  ends  of  the  layer  can  be 
written  as  follows: 
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M  =  0 
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where  M,  V  3-e  the  moment,  shear  and  a  is  the  average  value  of  the  slope 
-easu red  tmouab  the  thickness  given  bv 


A  suitable  elasticity  solution  that  represents  loading  on  the  upper 
surface  of  an  isotropic  elastic  layer  in  plane  strain  with  and  no  loading  on 
its  lower  surface  is  obtained  using  the  techniques  of  Sneddon  [A]  and  is  given 


as  follows: 
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where  3  =  ?h,  and  u,  v  are  the  shear  modulus  and  Poisson's  Ratio  of  toe  layer 
material,  respecti vely .  The  correspond! nq  equations  for  a  layer  in  piano 
stress  with  the  same  loading  conditions  are  obtained  by  replacing  v  with 
v/(Uv),  Fu  rt.hermore,  this  substitution  is  to  be  used  through  the  rest,  of 
this  paper  to  convert  plane  strain  expressions  to  corresponding  expressions  i 
plane  stress. 

It  is  seen  that  on  y  =  h  the  normal  and  shear  stresses  vanish 
automatically,  and  that  on  y  =  0,  the  shear  stresses  vanish.  The  normal 
stress  on  y  =  0  is  given  as 


=  /  [E<-(5)cos  Ux)  +  E^fOsinfCx)  ]d; 


12.16) 


We  let 
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Then  Eqn.  (2.16)  becomes 
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The  moment  and  shear  due  to  the  stresses  as  given  by  Eqns.  (2.11 )-(?.) 3)  are 
given  by 


Mix1  =  •'  vt  dy  =  -  f  !~Er  f  Tlcos (  ?:< )  +E ,  (^ )  si  n  ( Ex ) 1  — 
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We  note  that  Eqns.  (2.20),  (2.21)  satisfy  the  condition 


dM 

=  37  * 


(2.22) 


From  Eqns.  (2.10),  (2.15)  the  average  value  of  the  slope  is  found  to  be 
gi  ven  by 

OQ 

5E(x)  =  77  /  CEs(Osin(?x)  -  Ea(Oc os(5x)]d?.  (2.23) 

o 

For  the  beam  theory  solution,  the  displacement  is  taken  to  be  of  the  form 


ag  +  a^x  +  a2x2  +  a3X3 . 


(2.24) 


Assuming  the  hypotheses  of  Eu 1 er-3ernou 1 1 i  Beam  Theory,  the  following  results 
are  obtained: 
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Using  Eqns.  (2. 24 )-( 2.26) ,  the  moment,  shear  and  average  slope  are  calculated 
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where  0  =  uh3/6(l-u).  Thus,  the  solution  sought  will  be  a  superpos i 1 7  on  of 
the  elasticity  solution  given  by  Eqns,  (2.11 )-(2.16) ,  (2.201— (2.23)  and  the 
beam  theory  solution  given  by  Eqns.  (2.24)-(2.29 1 .  The  two  solutions  are 
matched  properly  wfien  the  boundary  conditions,  Eqns.  (2.1)-(2.9)  are  all 
satisfied.  This  is  achieved  by  solving  for  the  constants  a j ,  a 2  and  83  by 
superposing  the  two  solutions  and  applying  the  beam  theory  boundary 
conditions.  Thus,  Eqns.  (2.6),  (2.8)  and  (2.9)  become 


ME(4-i0)  +  ^B(^-^0)  =  0 

v£(t-i0)  +  =  0 

5g( -Zq )  +  (-Zfl )  =  0. 


(2.30) 

(2.31) 

(2.32) 
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Applying  Eqns.  (2. 20)-(2.23)  and  (2. 2d). (2.29 )  to  Eqns.  (2. 30). (2. 32)  yields 
the  following: 


.  <»  Zo(Z-Zo/2) 

aL  =  -  ^  | -  [Es(OsinC(Z-Z0)  -  EA(Ocos5(Z-z0)] 
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Zq 
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b  rt 


1_  "  S-( 3-2v)shS 
2y  1  S ( 3- s  n  S ) 
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(2,33) 
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- 2TT-  /  rEs(5)sin?fi-i0)  -  EA(5)cos5(A-Z0)l 


(2.34) 


aj  =  -ij  /  TEs(Osi  n5(i-iu)  -  Ea(Ocos5U-;Z.q  )  ]  . 

o 


(2.35) 


The  two  remaining  boundary  conditions  given  by  Eqns.  (2.5)  and  (2.7) 
still  need  to  be  satisfied.  First,  consider  Eqn.  (2.5): 


Uy (x.O)  +  uB(x,0)  =  A  -  x2/2R 


(2.36) 


Differentiating  with  respect  to  x: 


3uf  3uB 

rf(x-0)  *>r(x-0)  *  -i 


(2.37) 


Separating  Eqn.  (2.37)  into  a  symmetric  and  an  anti  symmetric  part  with  respect 
to  x  yields 


j~  (x,0)  |  +  2a2x  =  -  j 


(2.38) 


(x.O)  +  a^+  3a;jx2  =  0. 


(2.39) 


Consider  first  Eqn.  (2.38).  Substituting  for  tjr(x,0)  according  to  Eqn. 
(2.15),  combining  terms  in  E^(0  and  E^fO.  and  making  use  of  Eqns.  ^2.17) , 
(2.13)  yields 


/  *(t)  ' 


r  fh;f  3+sh3ch3  . 


°  32-Sd23 
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An  asymptotic  evaluation  of  the  kernels  in  Eqn.  (2.40)  shows  that  they 
are  all  convergent  at  the  lower  limit  of  integration.  However,  at  the  uppei 
l’.mit  the  first  term  in  the  first  kernel  is  divergent.  This  is  adjusted  by 
adding  and  subtracting  the  term 


h3 
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oo 

v(t)  I  sin(Cx)Jo(5t)d?dt 


0  0 

in  Eqn.  (2.40)  as  follows: 

-  r  I  +(t)  /  s1n($x)J0Ut)d5dt  +  /  «j,(t)  J  [!ll(g*.sh.gchg  +1 1  si nUx) 

oo  oo  S2-sh2s 
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C  00  x(Z-Iq) 

+  /  *(t)  /  [—  sin5(£-£0) - - -  cosSU-£0)]JiUt)d5dt  =  JL  (2.41) 

o  0  s2  5  R 

After  simplification  (i.e.  through  use  of  the  Weber-Schafheitlin  Integral 
[5]),  Eqn.  (41)  reduces  to 

h3  c  c  n 

■g—  'l'(x)  -  J  -4»(t)KA  (x.t )dt  -  J  ^ ( t ) K 2 ( x , t ) dt  =  -  ~  (2.42) 
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0  62-sh28  C2 
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<  2  (  X  ,  t  )  =  ~  t  x  . 


(2.43) 

(2.44) 


to  Eqn.  (2.39),  performing  similar  manipulations  v. elds 


“’eturni  ng 
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(2.45) 


An  asymptotic  analysis  reveals  that  all  the  terms  are  bounded  as  ?  +  0,  but  as 
r  +  *  the  first  term  of  the  first  kernel  is  divergent.  This  is  corrected  by 
adding  and  subtracting  the  term 


u3  C  ® 

g-  /  *(t)  /  cos(£x)Ji(5t)d£dt 


in  Eqn.  (2.45)  as  follows: 
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(2.46) 


simol i fi cation  Eqn.  (2.^6)  becomes 


(2.471 


•g-  4>(x )  +  /  i>(t)K3  (x  ,t  )dt  -  J  -?(t  )K4  (x  ,t  )dt  -  0 
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where 
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IT  2 
"  7  x 

(2.48) 

Kjx,t) 

=  *  hi  (S+sh3ch6  .  l)5xJl{5x)Ji(5t,d5. 

(2.49) 

o  32-sh23 

Sons.  (2.42)  and  (2.47)  are  the  two  coupled  integral  equations  tor  the  unknown 
auxiliary  functions  +(x),  <?(x).  These  equations  are  solved  numerically.  Once 
y(x),  $  ( x )  are  obtained,  all  necessary  physical  quantities  may  be  calculated. 
Stresses  may  be  calculated  using  Eqn.  (2.19): 
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(2.19) 


The  resultant  load  due  to  the  symmetric  stresses  is  obtained  as  follows: 


c  c 

P  =  -  J  t  dx  =  -IT  /  *(t)dt.  (2. SO) 

-c  yy  o 


The  resultant  moment  due  to  the  anti symmetri c  stresses  is  obtained  as  follows: 


_  c 

j  t<J»(t )dt . 
o 


(2.51 ) 


In  order  to  evaluate  the  deflection  under  the  indenter,  A,  the  constant  ar, 
must  be  determined.  Su peri moos: ng  Eqns.  (2.15)  and  (2. 241  yields 


uv(x,0)  =  ^  /  -■  — —  -  [Es(S)cosfCx)  +  E^(Osin(Sx)ld? 
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M  o  S(B2-sh2B) 
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(  2.52) 


Applying  Eqn.  (2.7)  yields 


a„=  -fill)  /i^!iSMrEc(5)Cos(5i„) 


E^( O si n ( ^  T-jr"  +  al20  "  a22$  +  a3^ 

(2.53) 
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where  a;,  a2,  83  are  given  by  Eqns.  (2.33)-(2.35) .  Substituting  Eqns.  (2.33) 
(2.35),  (2.53)  into  Eqn.  (2.52)  and  recalling  that  u  (0,0)  =  A  yields 
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(2.5A) 


An  asymptotic  analysis  reveals  that  both  kernels  are  bounded  as  5  0, 

but  the  first  two  terms  in  each  kernel  are  divergent  as  5  +  <n.  This  is 
adjusted  by  adding  and  subtracting  the  appropriate  terms  to  yield 

c  c 
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(2.55) 
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A  quantity  of  interest  in  this  problem,  as  well  as  in  the  subsequent 
study  of  the  dynamic  case,  is  the  rotation  of  the  beam  under  the  indenter. 
More  specifically,  the  rotation  of  the  beam  at  the  point  under  the  indenter 
about  which  the  moment  produced  by  the  antisymmetric  stresses  is  zero  (i.e. 
the  "eccentricity"  of  the  load)  is  of  special  interest.  This  point  is 
obtained  from  statics  simply  as 


x  =  -e  =  M/P. 


(2.58) 


Superimposing  Eqns.  (2.23)  and  (2.29)  and  making  use  of  Eqns.  (2.17),  (2.18) 
and  ( 2 . 33) -( 2 .35)  yields 
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(2.59) 


Note  that  the  boundary  condition  S(-£0)  =  0  is  automatically  satisfied  by 
Eqn.  (2.59).  An  asymptotic  analysis  of  the  kernels  in  Eqn.  (2.59)  reveals 
that  they  are  bounded  as  S  0,  but  that  as  ^  +  ”,  the  first  term  in  each  is 
divergent.  This  can  be  adjusted  as  before  to  become  the  following: 
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(2.60) 


Thus,  we  can  easily  evaluate  S  (x  =  -e). 


3.  PROBLEM  SOLUTION 


In  order  to  be  able  to  solve  for  the  auxiliary  functions  T f x ^  ,  i(x', 


-  ,n  n 


'2. a?' -(2. -A)  (2.A7M2.A3'  must  be  non-dimensi  onal  i  zed . 


Th is  is 


accomplished  by  the  use  of  the  following  non-dimensi onal  parameters: 
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y  =  x/c 

(3.1e,f) 

Once  the  non-dimensional i zed  auxiliary  functions  ¥fy),  $(y)  are  obtained, 

Eqns .  (2.191,  (2.501,  (2.511,  ( 2.55 )-( 2.571  and  (’.501  can  be  used  to 
calculate  non-dimensi onal i zed  stresses,  loads,  moments,  deflections  and 
rotations,  respectively.  These  non-dimensional  quantities  may  be  transformed 
back  to  real  quantities  via  the  following  relationships: 


T  =  P  T  /C 

yy  yy 

b  =  h2A/R 

(3.4),  (3.71 

p  =  h2uP/d-v)R 

5  =  h  §/R. 

(3.5),  (3.8) 

M  =  h3uM/(l-v)R 

(3.6) 

In  order  to  assess  the  accuracy  of  the  results  for  displacement  and 
rotation  under  the  indenter,  the  following  beam  theory  solutions  are  used: 

ab  =  3T)  P£^  +  In  n*9) 

"B  =  Til  Pf£o-e2)  +  ^  MUo-e)  n.in) 

It  shcjid  be  mated  that ,  although  the  formulation  of  the  equations  in 
t'nis  paper  is  sucn  that  tuey  can  he  adapted  to  plane  stress  or  mla^e  st-ain, 


the  use  of  beam  theory  equations  in  comparison  studies  dictates  that 
conditions  of  plane  stress  should  be  assumed  in  the  evaluation  of  Fans.  (?.55 
and  (2.60).  This  assumption  would  be  consistent  with  the  use  of  a  thin  hpam 
in  an  experimental  verification  of  the  results  presented  here.  For  the 
numerical  evaluation  of  Eqns.  (2.55)  and  (2.60)  under  conditions  of  plane 
stress,  Poisson's  Ratio  is  taken  to  be  equal  to  0.35  (this  corresponds  to  a 
Poisson's  Ratio  of  0.26  if  Eqns  (2.55),  (2.60)  were  to  be  evaluated  under 
conditions  of  plane  strain). 


i 


4.  OBSERVATIONS  AND  C0NCULSI0NS 


Solutions  to  the  problem  were  obtained  for  6  =  0.2,  0.5  and  1.0,  y  = 

5.0,  10.0,  15.0  and  20.0,  and  for  each  y,  a  =  0.25y,  0.5y,  0.75y.  It  was 
found  that  the  stress  distributions  for  the  different  values  of  1/ h  remained 
vi-tually  identical.  This  is  illustrated  in  Tables  1  and  2,  which  compare 
peak  symmetric  and  total  stresses  for  different  values  of  c/h  and  l o/h.  A 
typical  symmetric  stress  distribution  is  shown  in  Fig.  2;  that  of  the  total 
stress  under  the  indenter  is  shown  in  Fig.  3.  As  can  be  seen  in  Fig.  2,  for 
small  values  of  c/h  the  Hertzian  distribution  approximates  the  stress  under 
the  indenter  quite  well.  As  c/h  increases,  the  distribution  changes 
si gni ficant ly .  In  Fig.  3  we  see  that  as  c/h  increases,  the  location  of  peak 
st-ess  shifts  to  the  left  and  the  distribution  becomes  less  and  less 
parabolic.  This  is  because  for  small  values  of  c/h,  the  indenter*  acts  like  a 
coin*,  load.  He^ce,  for  small  c/h,  we  have  a  small  region  of  concentrated 

,  as  c/h  increases,  the  shaoe  of  the  indenter  becomes  more  ana 


st, -ess. 


-owever 


Tables  3  and  4  show  a  comparison  between  the  elast'c'ty  solution 


developed  here  and  the  classical  beam  theory  solution.  The  elasticity 
solution  is  seen  to  agree  well  with  the  beam  theory  solution.  However,  as  r/n 
approaches  i o/h  the  two  solutions  ditfer  to  a  much  1a"ger  extent.  This 
behavior  is  due  to  the  fact  that  the  elasticity  solution  is  valid  only  for 
c/h  <<  £o/h. 

It  is  ’important  to  note  that  the  rotations  in  Table  4  have  no  meaning 
unless  they  are  converter  to  real  values.  Since  the  elasticity  and  beam 
theory  solutions  were  derived  using  small  angle  approximations,  it  must  be 
ensured  that  =  <  5n,  where  3g  is  a  predetermi ned  angle  below  which  rotations 
are  conside~ed  “small."  Thus,  it  is  possible  to  place  a  lower  bound  on  R,  the 
radius  of  curvature  of  the  indenter: 


R/h  >  i/’Q.  (4.11 

However,  't  must  also  be  ensured  that  no  yielding  occurs  anywhe~e  in  the  beam, 

i . e . ,  shat  a  x  a  ,  where 
max  y 

Sna,  -  Mmaxh''2!»-  (4.?),  (4.3) 

and  a  is  the  yield  stress  of  the  material.  The  maximum  moment  occurs  at  the 
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fixed  end  of  tne  beam  and  is  given  by 
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R/h  >  6u( aP-M ) / ( 1  - vl  c 
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(4.5) 


If  =  20°  (0.34Q  rad),  =  55,000  ksi  and  p 
and  (4.5)  may  be  combined  into  the  following: 


11,510  ksi,  then  E q  n 5  .  (4.1) 


R/h  >  max{?.B653,  1.6a6(aP-M)} . 


(4.6) 


Judging  from  t'ne  values  of  3,  P  and  M,  it  75  evident  that  0.5653  >>  1.6P6(a 
for  each  case  considered.  Thus,  Fan.  (4.6)  reduces  to 


R/h  >  2.8653. 


(4.7) 


Under  these  criteria,  Table  5  shows  the  minimum  allowable  values  of  R/h 
to  ensure  small  deflections  and  no  yielding.  It  is  seen  that  for  large  c/h 
arc  for  long  beams  with  c/h  >  0.5,  the  indent er  is  practically  flat.  It  is, 
therefore,  of  some  interest  to  solve  the  problem  of  a  flat  indenter.  This 
solution  's  found  in  the  next  section. 


5.  CAST  II:  FLAT  INDENTER  -  PROBLEM  FORMULATION 


The  problem  to  be  solved  is  that  of  an  elastic  layer  of  thickness  h  and 
length  i  ■  n  den  ted  by  a  flat  punch  on  its  upper  surface  (Fig.  lb).  The 
'ourdary  ccndi t i ons  for  the  elasticity  problem  and  beam  tneorv  problem  vnose 
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tt  is  further  imposed  that  the  stresses  be  non-singular  at  x  =  c: 


Tyy(c,0)|  < 


(5.10) 


This  condition  ensures  a  smooth  deflection  of  the  beam  for  x  >  0  and  results 
in  a  “receding  contact"  [6];  here,  as  the  load,  increases,  the  contact  length 
decreases  with  load. 

Ths  suitable  elasticity  solution  for  this  problem  is  given  by  Eqns. 
(2.11)  through  (2.15).  Using  this  formulation,  the  normal  stress  is  given  by 


t  (x,0)  =  J  C  E<-  (  s)cos  5x  +  E,  (5)sinsxlds. 


yy 


(5.11) 


Tor  a  flat  pinch,  to  have  singularities  at  x  =  i  c,  we  let 


'-)  -  4,i,  (;c)  -  /  st)d- 


•  ■>] 
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c 

EaU)  =  BJ:Uc)  +  /  $(t)Ji(St)d5. 

0 


(5.131 


Substituting  into  Eqn.  (5.11)  yields 


A  +  i  B 

-  H(c-x) 

/c2-x2 


x  / 1 TT2  x 


The  first  term  in  Eqn.  (5.14)  is  singular  at  x  = 
vanish  at  x  =  c,  it  must  be  true  that  A  +  R  =  0. 


<p(t  )dt 


(5.14) 


tA2-x2 

±  c.  For  the  singularity 
Then  Fqn.  (5.14)  becomes 


to 


t  (x,0)  =  -  —  — H(c-x) 
yy 


+  j  4>(t)dt  +  x  j  4»  ( t  )dt  . 
x  A2-x2  x  t/t2-x2 


(5.15) 


The  moment,  shear,  and  average  slope  are  given  by  Eqns.  (2.20),  (2.21)  and 
(2.23),  respectively.  As  before,  the  beam  theory  solution  is  taken  to  be 


a  o  +  aix  +  a  2  x ' 


a3x~ 


(5.16) 


resulting  in  expressions  for  moment,  shear  and  average  slope  given  by  Eqns. 
(2.27),  (2.28)  and  (2.29),  respectively.  The  constants  a  ,  a^,  a^  are  solved 
for  as  before,  with  identical  outcome.  Eqn.  (5.5)  is  treated  in  a  similar 
manner.  Differentiating  with  respect  to  x  and  separating  into  symmetric  and 
anti  symmetric  components  yields 


( x  ,0) 

( x  ,0 ) 


+  2a2X  =  0 
A 


+  a  j_  +  3a3x2 
S 


0. 


(5.17) 


(5.18) 


Tonsidering  each  equation  ore  at  a  time  and  following  the  same 


'1'  J  !“* 
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as  before  results  ■  n  the  following  equations  for  4,  3 f  y(t)  and  t>(t): 


AK, (x,c) 


c  c 

+  ;  y(t)K; (x,t'dt  +  3K2(x,c)  +  /  5 ( t ) K2 f x , 1 1  dt 
0  0 


n  (5. IQ) 


AK  3  ( x ,  c )  +  /  :Ht  )K  3(x,t  )dt  -  AKjx.c! 


+  ^>(x) 


J  t)K„(x,t)dt  =  n  (5 .20) 


where 


oo  o 

'Ki(x,n)  =  J  {"7" — v - — — —  +  l},xJo(5xl  +  —  cosC ( o  1  ;0 o ( )d 5  +  i  x(4-1q) 

o  a2-s  h2a  ?2  * 

(5.?n 

*  J  nx  (5.22) 


K3 f x ,n)  =  -J  x2 


(5.23) 


<•+  ( x  ,n) 


f  h3  ' £+sh3ch3  + 

i  °  V-sh2£ 


lj£xJiUx)Ji(Cn)d£ 


(5.24) 


and  r,  =  c  or  t.  Eqns.  (5.19),  (5.20)  are  two  equations  in  four  unknowns.  A 
third  eauation  was  obtained  earlier  by  enforcing  non-singularity  of  normal 
stress  at  x  =  c: 


A  +  3  =  0.  (5.25) 

The  fourth  is  found  by  enforcing  the  boundary  condition  given  by  cqn.  (5,5). 
After  so! v' rg  for  the  constant  ao  as  before  (with  identical  results),  cqns. 
(5.5),  (2.15),  (2.24),  (5.12),  (5.13)  and  (5.25)  may  be  used  in  a  manner 
ana! conu 3  to  that  so!  •/:  no  for  t~°  d*  sol  acement  under  the  cyl  i  ndni -3 1 
!  n d n t  e “  to  s o ’  v e  ror  t h -■*  c n -  - 1 0  - 1  ; 


(5.26^ 


^  -  j  Y(t)K5(t)dt  -  /  p(t)K6(t)dt 
0  0 
KiTcT  -  K 5 1  c  ) 


£7  1(B^h£t 

o  e2-sh2e  4 


sinC^o  3z§  cosCf l-lo ) 

£°  FsW  '  ‘  ^ - ~2 - }  Jo  f^nldC 


,  TT  3-2v 
1  ~ uTT 


lQ 


~  Y2]J  3i~z g ) 


(5.27) 


3+sh6chS 

02-sh2S 


1) 


si nS* o 
l 


cosCio 


1-v  n  trri 

U  j -  “  BU 

Zq+1 ijj-n2 


(5.28) 


or  t.  Making  use  of  Eqns.  (5.25),  (5.26)  in  Eqns.  (5.19),  (5.20) 


C  <i  (x  ,C)-I<2  (X  ,C) 

»(x’  -  l  *ftu<‘(x"t)  +  -<7(c)  -  k5(c)  Mtlldt 

c  <i(x,c)-K2(x,c) 

-  /  *(t)!K2(x,t)  +-r-T-1_T-T_r  K6(t);dt 

-<1  (x  ,c)+K2(x ,c) 
=  A  <6(c)  -  K5(c) 


c  <3(x,c)+K4(x,c) 

?(xl  *  /  t(tHK3(x,t)  *  -^ic)  .  k5)c)  Mt»dt 


c 

-  /  ?(t){Ku(x,t) 


u 


<3  (x ,C )  +  <4 (x  ,C  ) 

"K6TcT  -  k  5 ( C  )  K6(t)}dt 


(5.29) 


<  3  (  X  ,  O  +<  ^  (  X  ,  c ) 


Thus,  Eqns.  (5.?9),  (5.30)  are  those  used  to  solve  for  the  unknown  auxiliary 
functions  y(x),  ;(x).  These  equations  are  solved  numerically.  Once 
y(x),  ?(x)  are  obtained,  all  necessary  physical  quantities  may  then  be 
cal cul ated . 

Stresses  may  be  calculated  using  Eqn.  (5.15): 

T  (x  ,0)  =  -  ~  -  H(c-x)  +  /  +  X  /  -  (5.31) 

/ :2-X2  X  /t2-X2  X  t/t2-x2 

where  3  is  given  by  Eqn.  (5.36). 

The  resultant  load  due  to  the  symmetric  stresses  is  found  to  be  given  by 

c 

p  =  rr8  -  7T  /  *(t)dt.  (5.3?) 

0 

The  resultant  moment  due  to  the  antisymmetric  stresses  is  found  to  be  given  by 

M  =  -  £  C8  -  j  j  1 4> ( t  )dt .  (5.33) 

o 

The  rotation  of  the  beam  at  any  point  x  >  c  is  found  to  be  given  by 

-  ?  c  c 

S(x)  =  B[<3(x,c)-'<7(x,c)]  -  B  -J  +  /  'Kt)K?(x,t)dt  +  J  *(t)<8(x,t)dt 

o  n 

(5.34) 

where 

cosy (4-t0 ) 

(x+Zn) - 1  Jn^nldr 

,  j  -j 


v  .(x  =  c. 


1  , 
J  j 


h3/sinyx+siny£^ 


-  sr 


-  ?  4  - 


(5.35) 


-  Jry  [io(^^o/2)+x(l-i0)-x2/3]  +  *  li2v 

.  «  3  COSSx-COS^o 

'<Q(x,n)  =  iy  j  -g— iv — I- sFT5 - jJ  i(  ?n)d^  -  (io+x^r'  (5.36) 

o 

and  n  =  c  or  t. 

6.  PROBLEM  SOLUTION 

Following  a  scheme  similar  to  that  of  Sec.  3,  we  define 

a!  =  *L/h  *(x)  =  ^  f(y)  (6.1),  (6.2) 

h5 

*(x)  =  *(y)  (6.3) 

h5 

Non-dimensi onal i zed  stresses,  loads,  moments  and  rotations  are  calculated  and 
may  be  transformed  back  to  real  quantities  through 


T 

=  Pt  /c 

(6.4) 

yy 

yy 

p  = 

uAP/( 1 - v) 

(6.5) 

M  = 

uAhM/(l-v) 

(6.6) 

W 

d  = 

A'/h. 

(6.7) 

In  order  to  assess  the  accuracy  of  the  solution  obtained,  two  tests  are 
pe-formed.  First,  the  rotation  of  the  beam  just  beyond  x  =  c  is  compared  to 
t'^at  given  by  the  standard  beam  theory  solution: 
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where 


e  =  c  .  Second,  from  the  nature  of  a  receding  contact  problem,  as  the 
contact  lengtn  decreases  to  zero,  the  load  increases  to  a  certain  limit 
value.  Because  the  moment.  goes  to  zero  as  the  contact  length  goes  to  zero, 
the  beam  theory  solution  for  displacement  tells  us  that  the  load  must 
approach  the  beam  theory  load  Pg^  for  a  cantilever  beam  problem: 


o 

'  3T 


30A./4C  . 


(6.9) 


The  "limit  lead  lim  P 
c/h+0 


compared  to  the  value  Pgr  for  each  case 


consi dered . 

Once  again,  the  use  of  beam  theory  equations  in  subsequent  comparison 
studies  dictates  that  conditions  of  plane  stress  should  be  assumed  in  the 
solution  of  Eqns.  (5.29)  and  (5.30).  As  before,  Poisson's  ratio  is  taken 
be  equal  to  0.35  for  the  numerical  solution  of  these  equations. 


to 


7.  OBSERVATIONS  AND  CONCLUSIONS 

Solutions  to  the  problem  were  obtained  for  5  *  0.25,  0.5  and  1.0,  y  = 
10.0  and  20.0,  and  for  each  y,  i  =  0.25y,  0.375y,  0.5y  and  0.75y.  It  was 
found  that  the  stress  distributions  for  the  different  values  of  1/ h  were 
virtually  identical.  A  typical  distribution  of  the  total  stress  under  the 
"reenter,  normal ized  by  a  factor  of  /c 2 - x 7 ,  is  shown  in  Fig.  5.  An 
examination  of  the  plots  of  the  total  stress  under  the  indente-1  reveals  that 
it  is  indeed  singular  at  x  =  -c,  but  zero  at  x  =  -c.  Furthermore,  t~e  stress 

^  *  •"  ■j  ^  •  ’  *  :t  ^  rj  *  rj  ~  #  5  f  C  T  C  /^*  =  i  •  0  ]  3  *■  S  r  n  "omp:  £  ~i  <b  3  ^  ^  i  S 
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results  from  a  more  pronounced  effect  of  the  antisymmetric  component  in  the 
solution  for  cases  where  c/h  is  large. 

Tables  6  and  3  show  a  comparison  between  the  solution  developed  here  and 
the  classical  beam  theory  solution.  As  can  he  seen,  the  elasticity  solution 
agrees  quite  well  with  the  beam  theory  solution.  Again,  it  is  noted  that  the 
rotations  in  Table  6  have  no  meaning  unless  they  are  converted  to  real 
values.  Following  the  same  procedure  as  before,  an  upper  bound  can  be  placed 
on  the  ratio  of  displacement  under  the  indenter  to  beam  thickness: 

A/h  <  min  {o0/3  ,  ( l-v)a  /fiy(aP-M) }  (7.1) 

Assuming  the  same  values  of  3o,  and  u,  Eqn.  (7.1)  reduces  to 

A/h  <  0. 349/3.  (7.2) 

Under  these  criteria,  Table  7  shows  the  maximum  allowable  values  of  A/h  to 
ensure  small  rotations  and  no  yielding. 

The  phenomenon  of  receding  contact  in  this  problem  is  borne  out  by  Fig. 
7.  As  c/h  goes  to  zero,  the  non-dimensi onal i zed  load  parameter  approaches  a 
limit  value  greater  than  zero.  In  Table  3,  the  limit  load  of  the  elasticity 
solution  is  seen  to  agree  quite  well  with  the  limit  load  of  the  beam  theory 
solution.  Furthermore,  the  limit  load  is  highest  in  those  cases  where  the 
indenter  is  close  to  the  wall. 

In  Fig.  7  the  relationships  between  the  loads  and  contact  widths  for  the 
va-ious  cases  considered  appear  to  violate  the  condition  imposed  on  receding 
contact  o-ob’ -ms  to  e°s.:re  the  e<i  s*- ^nc-3  and  uni oueness  of  fhei~  solutions 
y- ' ,  s i  '"te,  tne  contact  1  e-ctos  *or  a  parti  cu 1  ar  beam  are  not  independent  of 


the  applied  loads  on  that  beam.  The  contact  length  in  a  receding  contact 
problem  will  be  independent  of  the  level  of  loading  only  i *  the  ra^.io  of  tho 
resultant  moment  to  the  apolied  load  is  constant  '"or  all  levels  of  loading. 
In  the  problems  solved  through  the  course  of  this  stu^v,  tHe  ratio  of  moment 
to  load  changes  as  the  applied  loading  changes;  thus,  the  contact  length  1 
change  as  well,  the  imposed  conditions  are  not  violated,  and  the  problems 
studied  are  indeed  receding  contact  problems. 
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Table  3:  Displacement  Comparison  (Case  I) 
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Table-  6:  Rotation  Comparison  (Case  17 
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An  investigation  i s  made  on  the  low  velocity  impact  on  a  ci  rcu  1  ar  plate 
by  a  r-'g'd  indenter.  The  plate  response  to  impact  shows  some  of  bat  not  all 
of  tne  similar  problem  for  a  beam.  The  local  contact  stresses  were  developed 
from  an  e’asticity  theory,  wnile  tne  global  stresses  were  obtained  from  plate 
theory . 


Tne  problem  of  the  low  velocity  indentation  of  an  isotrop''c  circular 
elastic  plate  by  a  rigid  purer  has  been  'investigated.  This  study  uses  the 
assumption  that  an  elasticity  tneory  may  be  used  for  calculating  the  local 
effects  and  a  Bernoulli:  -  Euler  plate  theory  can  be  used  for  calculating  the 
overall  dynamic  response. 

In  tn-’s  problem  a  spherical  indenter  with  a  very  large  radius  of 
curvature  strikes  the  plate  with  a  low  impact  velocity.  Since  both  the 
pressure  load  and  contact  area  are  unknown  and  depend  upon  time,  the  problem 
is  formulated  and  reduced  to  solution  of  a  Volterra  integral  eguat'en. 

”wo  types  of  solutions  are  superposed.  One  is  the  static  solution 
analogous  to  the  one  developed  by  Keen  and  Miller  gi];  the  other  is  a  standar 
dynamic  plate  solution  whicn  can  be  obtained  from  Sernoui 1 1 i  -  Euler  plate 
analysis  !see  e.g.  Graff  [2]). 


re  pres  >  i  re  d*'  stri  out 4  on  is  assumed  to 


?w'se  cons; ant 


o  r  a'  i  n  :  r  :*n  static  and  dvnem:  ■;  so!jt'-rs,  the  prws  Sire  d’  Str:  but  m 


wi:  ;i;a  trat  was  ce/e  1  ope  ; 


a  n  1  ■  > , 


tie  time  wise  pressure  distribution  's  calculated. 


An  example  for  the  case  o 


plane  strains  low  velocity  has  been  given  in  an  earlier  paper  t>y  Keer  and  Lee 

[A]. 


t 


PROBLEM  FORMULATION 

The  formulation  uses  the  assumption  that  el astodynami c  effects  are 
negligible  wnen  compared  with  the  effects  due  to  the  static  deformation, 
nowever,  dynamic  effects  are  incorporated  into  the  plate  bending  solution. 
Tne  solution  of  this  problem  is  obtained  by  superposition  of  the  static 
solution  with  a  dynamic  plate  solution. 

The  boundary  condition  for  an  elastically  supported  layer  loaded  by  a 
smooth  rigid  indenter  are  as  follows: 


tzz  (r,h)  =  0  0  <  r  <  °°  fl) 

(r,h)  =  0  0  <  r  <  °°  (2) 

-zr(r,o)=0  0  <  r  <  ™  (3) 

'zz  (n,o)  =0  c  <  r  <  »  (4) 

( r ,o )  =  -p(r)  0  <  r  <  C  (5) 


where  c  Js  tne  radius  of  the  contact  area. 

'ne  boundary  conditions  at  the  ends  of  the  plate  are: 


uz=0  r  =  a  (6) 

M  -  <"  =  0  r  =  a  (7) 

wre'-e  a  '■  3  t he  radius  of  tne  plate,  M  is  the  moment,  <  is  spring  constant, 
:  '  s  tn®  value  0f  tne  slope  averaged  through  the  thickness. 


he  elastic:  tv  solution  that  represents  loading  on  the  upper  surface  of 


e'ast't  1 3yer  is  given  in  the  following  for n  by  jse  of  integral  transforms 

2uu „  =  : (l-2v)(3  +  Sh3ch3l  -  ^zsh 2  3 ' sh ( z ) 

-  _(l-2vlsh23  +  32  -  £z(3  +  shsch3)ch(;z) } 

*  Ji(;r)d^  ( 9 ) 

2uu  =  -  ”  E4~~~  i  2(l-v)sh23  -  3 2  +  -;z 

Z  0  o 

(3  +  sh3ch3;sh(-;z)  -  -2(l-v)(8  +  shschsl 

-  Czsh23]ch(rz) } J Q ( r r ) d C  (10) 

t  =  zE-X--  l[sz  -  8  -  sh  3ch  3 ] sh  Uz) 
rr  o  a  k  L 

f  '32  +  sh23  -  5Z(3  -r  Sh3ch3)]ch(-z)  jJ0(cr)d' 

-if"  iiiii  { [ 3 2  +  (1  -  2v)sh2s  -  Zz 
[3  +  sh  6ch  3) ]ch (£z)  -  frZsh23 

-  ( l-2v)  (3  +  shSche) Ishf'Z) }Ji (Cr)d^  (11) 

T  =  !*“  r  E ’’  { i’ 32  -  Zz  (3  +  Sh3ch3)  nsh(rz) 

r  £  o  A  1  L  J 


zs h 2  3  ■  ch  ( z  z )  ’*  J  i  ( -  r ) d ^ 
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,  ,  j  o  modulus  arc  •'o:  uon  s  nai :o, 


'.,*er-  ►  ■,  "  ,  ^v-'j  . ' ::  s'ne  an  :  cos' ne,  respect' ve !  y. 

:  s-v^  st '•oss  /T°;sn  at  z  =  n ,  but  at  z  =  0,  the  shear 


?ss  : s  zoro ,  t-e  oor-a ;  s"'o< ;  's  r /on  as 


,  ,  r  ,  o 


r  •.,  •  n  ' r 


Taking  the  quadratic  di $pl accent  {constant  moment  solution)  as  in  the 
ate  solution  by  Keen  and  Vilier  ’1], 


u  3  Ar  *■  t  B 


wnere 


-  \(L)  Eh3 

23(7  *  v)  ’  J  T2T1"^T  11 

*o  determine  tne  constant  A,  we  first  derive  the  moment  and  average  slope 
ecosentati  on  due  to  the  layer  solution  in  such  a  manner  that  the  end 


From  equations  (10),  (lb),  (16),  120',  and  (221 


1  -c  ■co  i  h3  3  +  shdcns 

2D  o  K  oo  l3  FstvDT 

«  (J0(sa)  -  J„(cr;)  -  [ispi 


1  ,  ,  ,  0  ,  32ch3  -  sh 2 3,  s  ,  , „ ,  N 

IF  (1  +  2v  (-TTW))J-(;a) 


+ 


h2a 

I  2(1  -  v ) a 


(3  -  2v ) sh  3 
3  -  sn  3 


)Ji 


5a) 1  Id: 


(23) 


As  ;  +  0,  an  asymptotic  analysis  shows  the  terms  in  equation  (23)  to  be 
bounded.  However,  as  E,  +  ®,  the  terms  are  divergent.  By  adding  and 
subtracti ng 


^  oc  r'p(r')tir'  £*  (JQ(ca)  -  JQ(5r))d5  (24) 

„  (3-2v)  a2  -  r2  h2a  rc  .  /„*n 

x  "”2D~  r~  VI  12(1  '-“via  o  r  p(  } 


*  fir  J0(5r")J1 (ca)d^dr^ 


(25) 


tne  difficulty  can  be  removed. 


DYNAMIC  PLATE  THEORY  SOLUTION 

A  solution  can  be  obtained  by  using  a  normal  mode  expansion  together  wit 
3  Laplace  transform.  The  equation  governing  the  normal  mode  expansion  is 
; :  / e r  as 


r 


Since  both  V : ,  '<:  nave  a  singularity  at  r  =  0,  they  must  be  set  to  zero 
Tn us,  the  solution  is  reduced  to 


Yn(r)  =  i A J  o ( 3  n  r )  +  C I o ( 3n  r ) }  (31) 


For  tne  clamped  case 


av  (a  1 
n  ' 


y  iai  =  ^  '  = 

0 

(32) 

n  3r 

Substituting  into  equation 

(31),  we  get  the  frequency 

equati on 

I  a ( i  d ) Jo ( 3„a )  - 

3  3  (  =na  )  I  3  (  3n<3  )  =  0 

(33) 

apd  the  normal  mode 

Jq(3  a) 

•  TyrTTy 

Yn(r)  =  iJo(3nr) 

(34) 

For  the  simply  support 

ed  case,  the  boundary  cond: 

Itions  are 

v  (a)  =  0  at 

r  =  a 

(35) 

n 

M  =  0  at 

r  =  a 

(36) 

oar/ 


;tions  35)  and  ( 36)  tne 


u  ency 


and 


'O^m  a 


-  1 1  ( 3  a )  J :  (  3  a )  =  0 
n  * '  n 


Y  =  J .t (  i  r)  - 


Jg ( 3  a ) 


);1(3  r) 


To  derive  the  dynamically  loaded  solution,  we  take  the  Laplace  transform 
of  the  governing  equation, 


D(ypr  +  7  fr)2Y(r)  +  "h  ~  =  P^’^ 


wnere 


rne  Laplace  transform  of  (39)  is 


32(11,  +  1  l_)2?(r,s)  +  s2 7 ( r  , s ’ 
Tr2  r  Ar  ’  ' 


p(r,s)/ph 


and  tne  transform,  7(r,s),  p(r,s)  may  De  expanded  in  a  series  of  normal 


modes  V  (r)  as 
n  x 


7(r,S)  *  i,  3_(s)Y„(r)  , 


ir7 11  + 


J  1  (  £3  ) 


h2a  ,3  -  (3  -  2v)sh3x  , 

+  T7TT  -  vIF  ( - 3“-_sh3 - }  Jl(’a)jtdc>dr 


-  1  1 
n  =  l  pF  UP 


Y  (r) 
n 


+  j  n  (  d  S 


x  ;'a  r'Y  (r')  p(u,0  sin  G62  (t  -  £)d*dr'  (51) 

5  n  o  n 


In  the  contact  reqicn  we  have 


(  r  ,  t  )  -  6  ( t  )  -  yrr  0  <  r  <  c  ( t ) 


wnere  ~:(t)  is  the  relative  approach  and  R  is  the  radius  of  curvature  of  the 
striker. 

From  Newton's  equation  of  motion 


5(t)  =  Vt  -  -4^2-  jt  (t  -  i)  fC^J  rp(r,r)drd5  (53) 

pna  o  'o 


wne-e  V  ■ s  the  velocity  of  the  striker  as  it  hits  the  plate  and  M  is  the  rati o 
r  tne  'n-:ss  of  the  plate  to  that  of  the  stri<er.  The  method  to  determine  the 


JN  AND  RESULTS 


r*  - 


ViRICAL  SOL 

In  solving  tne  Integral  equation  (52),  some  techniques  have  been  used  to 
ta«.e  care  of  the  singularity  and  speed  of  convergence.  Th;rd  order 
polynomials  are  used  to  approximate  the  unknown  function,  since  increasing  the 
order  of  polynomia’  to  higher  order  does  not  mean  that  better  accuracy  will  be 
achieved.  Filon's  formula  is  used  to  improve  the  accuracy  of  the 
tr-1  gnometri  cal  parts  of  the  kernal  .  The  overestimated  contact  area  and 
rod'Ified  starting  value  of  the  Volterra  integral  equation  help  to  ensure  the 
acc^acy  of  the  numerical  results. 

Calculations  are  based  on  a  variety  of  geometrical  and  material 
parameters.  They  are  the  following: 


Radius  of  Strixer: 
Radius  of  Dlate: 
Elasticity  Constant: 
rhic<ness  of  Plate: 
^oisson's  Ratio: 
Density  of  Plate: 


40  cm 
1  0  cm 
200  Gpa 
1  cm 
0.25 

7.8335  g/cm3 


figures  2-2  show  the  contact  area  and  contact  force  for  cases  of 
different  velocities  and  mass  ratios.  As  tne  mass  ratio  is  increased  the 
oe-iod  o‘  impact  becomes  shorter;  however,  changing  the  initial  velocity  does 
"a t  appear  to  a* feet  the  per: od  of  impact.  8y  comparing  with  Keer  and  Lee 
f’e  „er*od  o"'  impact  for  the  plate  is  shorter  than  that  of  the  beam, 
s  •  i  •:  ■nside'-ed  to  be  due  to  tne  stror>ger  boundary  effect  in  the  three- 


Tne  pressure  distribution  is  shown  in  Fig.  4-7,  which 


prominent  os  'o  [4], 
corresponds  to  several  values  of  the  mass  ratio.  The  non-rlertzi  an  solution 
oeg'ns  to  appear  in  Fig.  4,  and  the  wrapping  effect,  Keer  and  Miller  [1], 
becomes  more  important.  The  pressure  distribution  begins  to  concentrate  more 
closely  to  the  edge  of  the  contact  than  the  center. 

Figures  8-10  show  the  pressure  distribution  for  different  velocities, 
it  is  seen  that  the  ratio  of  the  peak  forces  tends  to  be  proportional  to  the 
ratio  of  velocities.  Ths  shape  of  the  force  history  and  contact  area  are  very 
s'nilar  for  these  cases. 

Figure  11  shows  the  pressure  distribution  for  the  impact  of  a  clamped 
plate  by  a  striker  with  velocity  20  m/s  and  a  mass  ratio  1.  As  can  be  seen 
the  rate  of  increase  of  the  force  is  the  same  as  in  the  case  of  the  simple 
support  plate.  Non-Hertzi an  solution  occurs  during  a  portion  of  the  impact 
time.  Tne  pe-iod  of  impact  is  noticeably  shorter  since  tne  clamped  plate  is 
s  t  i  f  f  e  r . 

Figure  12  shows  the  force  history  and  contact  length  when  the  thickness 
of  tne  plate  is  0.125  cm.  As  the  plate  becomes  less  stiff,  we  observe  a 
double  scrike  phenomenon,  similar  to  that  noted  by  Keer  and  Lee  [4]  for  the 
related  beam  problem.  The  velocity  of  the  striker  fortties  case  has  been 
reduced  to  2  m;s  to  make  elastic  behavior  assumption  valued  as  the  thickness 
of  tne  plate  is  reduced.  Figures  13,  to  14  show  the  pressure  distribution  for 
tnis  reduced  tnicsness  case.  At  the  timewise  point  near  detachment,  the  plate 
ex-erienced  is  rather  high  stress. 


•JW 


de  a  -et  ry  .  However,  toe  dou  d  1  e  stri  <e  phenomenon  con  still  occur  if  tn«  pi  ate 
is  sjt'fi  ciently  sett.  Related  observations  to  our  results  are  calculated 
Goldsmith  .  For  M  -  0.75  and  ^e]  vel  y  thick  plates  t°e  relation  “^w-en 
force  and  tine  ;s  noted  to  be  nearly  constant  for  a  portion  of  the  t me . 
Parameters  such  as  these  may  prove  to  be  useful  for  the  design  of  an  impact 
test  in  wh'cn  the  force  is  kept  relatively  constant  during  a  prescribed  period 
)f  time. 
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Seometry  for  Circular  Dlate  Loaded  by  Pigid  Indenter 

Contact  Lenth  Vs.  Time  for  Simply  Supported  Plate  Subjected  to 
Impact  ( v  =  20  m/s) 

Force  History  for  Simply  Supported  Plate  Subjected  to  Impact 
(v  =  20  m/s} 

Pressure  Distribution  for  Simply  Supported  Plate  (M  =  0.75, 
v  =  20  m/s) 

Pressure  Di st ri but’ on  for  Simply  Supported  Plate  (M  =  1.0, 
v  =  20  m/s) 

Pressure  Distribution  for  Simply  Supported  Plate  (M  =  1.5, 
v  =  20  m/s) 

Pressure  Distribution  for  Simply  Supported  Plate  (M  =  2.0, 
v  =  20  m/s) 

Pressure  Distribution  for  Simply  Supported  Plate  (  M  =  1.0, 
v  =  15  m / s ) 

Force  History  for  Simply  Supported  Plate  Subjected  to  Impact 
(M*  1 .0) 

'wUntdcu  vd.  i  i  me  i  vji  icu  r  iaoc  juujcutcu 

Impact  (M  =  1 .0) 

Pressure  Di  stri  buti  on  for  Clamped  Plate  (v  =  20  m/s, 

•1  =  1  .0) 

Contact  Lengtn  vs.  Time  for  Clamped  Plate  Subjected  to 
Transverse  Impact  (v  =  20  m/s,  M  =  1.0) 

-orce  History  for  Clamped  plate  Subjected  to  Transverse 
Impact  (v  =  20  m/s,  m  =  1.0) 

Force  History  for  Simply  Supported  Plate  Subjected  to  Impact  ( 
2  m/s.  Thickness  of  Plate  =  0.125  cm) 

Contact  Lengtn  vs  Time  for  Simply  Supported  Plate  Subjected  t  j 
Impact  7  =  2  m/s.  Tn’CKness  of  plate  =  0.125  cm! 

-'ressure  Distribution  for  Simply  Supported  Plate  ( V  -  1.0,  «  - 
i/s.  Tn'Ickness  of  3late  =  0.125  cm) 

-  re ;  >  u  -e  Di  str-  out  ion  •  or  S’moly  Supoorted  Plate  '■'*  =  1.5.  v  - 


1.0  GPa 
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